global  phii vex vey dy1p5 dx1p5
%dndz=sbdel(deni);
%dphidz=sbdel(phi);  
vey(2:nx-1,2:ny-1,:)=dx1p5*(phii(3:nx,2:ny-1,:)-phii(1:nx-2,2:ny-1,:));
vex(1:nx-1,2:ny-1,:)=dy1p5*(phii(1:nx-1,3:ny,:)-phii(1:nx-1,1:ny-2,:)); %free- fix


vex=vex;
vey=-vey;
vey=sbcy(vey);




